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■ Abstract: Factor models have large potencial in the modeling of several natural and human phenomena. In this 

paper we consider a multivariate time series Y„, n > 1, rescaled through random factors T^, n > 1, extending 
some scale mixture models in the literature. We analyze its extremal behavior by deriving the maximum domain 
of attraction and the multivariate extremal index, which leads to new ways to construct multivariate extreme value 
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. distributions. The computation of the multivariate extremal index and the characterization of the tail dependence 

show the interesting property of these models that however much it is the dependence within and between factors 
T„, n > 1, the extremal index of the model is unit whenever Y„, ti > 1, presents cross-sectional and sequencial 
tail independence. We illustrate with examples of thinned multivariate time series and multivariate autoregressive 
processes with random coefficients. An application of these latter to financial data is presented at the end. 

. Keywords: multivariate extreme value theory, factor models, tail dependence 
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1 Introduction 

Factor models have been used in the modeling of data within hydrology (Nadarajah [551 [27] 2006/2009, 
Nadarajah and Masoom [55] 2008), storm insurance (Lescourret and Robert, [52] 2006), soil erosion in 
, crops (Todorovic and Gani [M] 1987, Alpuim and Athayde [5] 1990), reliability (Alpuim and Athayde [5] 

1990, Kotz et al. [19] 2000), economy (Arnold, [3] 1983) and finance (Ferreira and Canto e Castro, [13] 
^ ; 2010). 

cn 

, Let X„ = (Xni, . ■ . , Xnd)^ n > 1, be a d-variate sequence, such that X„j = YnjTnj, j = 1, 

' where 

m 

(a) Y = {{Yni, . . . ,Ynd)}n>i IS a stationary sequence such that, Ynj has a Pareto-type distribution 
Fy^ , J = 1, . . . , d, i.e., for each j — 1, . . . ,d, there exists a positive constant f3j for which 

><: FY,{x) = l-x-f'HY,{x), (1) 

- - -' with a slowly varying function, i.e., ly. (ax) / ly. (x) 1, as a; — >■ cxd, for all a > 0, 

(b) T = {{Tni, . ■ . ,Tnd)}n>i IS a stationary sequence, independent of Y, with support M'J_ and such 
that E{T^^) < oo, j = 1, . . . ,d. 

This work is concerned with the extremal behavior of the multivariate time series X„ , extending most 
of the factor models mentioned above. More precisely, we derive the max-domain of attraction (Section 
|5|) , calculate the multivariate extremal index (Section and characterize the tail dependence (Section 

u. 

The product YnjTnj can be seen as a random normalization of Ynj by Tnj, which is often required 
when modeling extremal behavior. For instance, if Ynj is the rate of an extreme event and Tnj its average 
cost, then YnjTnj can be interpreted has the total cost of the extreme event. Products of two independent 
random variables where one of them is regularly varying have been addressed from both theoretical and 
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applied points of view (Maulik et al. [H] 2002, Lescourret and Robert [H] 2006, Nadarajah [H] 2006 
and references therein). 

Our motivation to the probabiHstic study of extremes of multivariate sequences of products was 
originated from some particular models. Consider, for instance that T„j are Bernoulli distributed. Then 
X„ provides a model for multivariate data subjected to missing values. Extremes of univariate sequences 
with random missing values have been considered in Weissman and Cohen (|35). 1995) as a particular 
case of some mixture models. Additional results on extremes of incomplete samples can be found in 
Mladenovic and Piterbarg ([55], 2006) and Zhongquan and Wang ([37], 2012). 

Li ([121 J 2009) analyzed the tail dependence of the scale mixture X„ when Y„ = ■ ■ ■ , Ynd) has 
multivariate extreme value distribution with standard Frechet margins and Tnj = T„ , j — 1 , . . . , d. Here 
we consider scale mixtures of multivariate sequences which are very flexible models for data exhibiting tail 
dependence and asymptotic tail independence such as, respectively, ARMAX and pARMAX sequences 
(Ferreira and Ferreira |15|). We give particular emphasis to a model in which /3j = a/"/j, a,7j > 0, 
j ~ 1, . . . , d, generalizing the results of Lescourret and Robert ([22, 2006) (Section[5|). An application to 
financial data will be provided at the end (Section |6|). 

2 Preliminary results and max-domain of attraction 

We start with some properties of {X„}„ > 1, that will be used along the paper. 

Proposition 2.1 For each j ~ I, . . . , d, {Xnj}n>i *s o, stationary sequence having Pareto-type distribu- 
tion. 

Proof First, observe that 

P(y„ Tr, >X) / {^l^)dPT„, iz) 

lim ^ — — 7 — = lim — — — = r,- , 

x^aa P{Ynj > Xj x-^oo x PHYj[x) 

where the last step is due to the dominated convergence theorem and by using the Potter bounds of 
regularly varying functions (Bingham et al.. [5] 1987; Theorem 1.5.6.). Therefore, for large x, 

\-Fx,{x)^P{X,,j>x)^x-l'HY,{x)r,{l + o{l)) := x^^Hx^x), (2) 

where it is immediately seen that Ixj is a slowly varying function. □ 

In the sequel we denote Uxj{x) and C/y^ (x) the quantile functions, F^^{1 — l/x) and ^^"'^(1 — 1/a^), 
respectively. 

Given ([1]), we can state 

Uy,{x)=x''PHu^Ax), (3) 
where lu^- is a slowly varying function, and by Proposition 12. 11 we can also write 

Ux,{x)=x^/^Hu^^ (x), 

where lu^, is a slowly varying function. By the Bruyn conjugate concept (Beirlant et al. [J 2004, 
Proposition 2.5), we have that, for large x, 

Ux, (x) = x^/^^ (x^/^^ lu.^ ix))il + o(l)). 

By (12) and we have for large x 

Ux, (x) = x^/^Hu^^ {x)r]'^^ (1 + oil)) = Uy^ (r, x) (1 + o(l)). (4) 
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Proposition 2.2 The upper tail copula function ofK. is given by 



Ax(xi,...,xd) = £; Ay 



, . . . , 



with (xi, Xd) G ~ [0, oo]''\{(cxD, . . . , oo)}, where T = (Ti, . . . , T^) is a random vector distributed as 
T„ — (r„i, . . . ,Tnd) and provided that the upper tail copula function o/ Y„ exists, i.e., the limit 

AY{xi,...,Xd) = ^lim <P l^fl {Yij > Uy, {t/xj)]^ (5) 

is finite. 

Proof The upper tail copula function of X is defined by 

K^{xi,...,Xd) = Imi^tP l^pl {Xij > Ux, {t/xj)] 
The result follows immediately by applying Q and the dominated convergence theorem, since 
}}^tP\T\{X,,>Ux,{tlx,)) 



rf^(Tii,...,Ti^)(^l,---,Z<i) 



(6) 



Remark 2.3 If Pj = j3 and T„j = T , j = I, . . . ,d, then Ax = Ay, corresponding to the case considered 
in Li (J^^j, 2009). If Ay{xi, Xd) = X^iLi^^i; ''^^ '^'■^o obtain Ax = Ay, for any choice o/T. In the 
last section we give particular attention to the case Ay(xi, ...,Xd) = AiLi ^j- 

Normalized levels u^Jj'^ of Xj are such that n(l - Fx^iul^f)) Tj > 0, as n — ^ oo, I.e., u„ = 
Uxj (n/Tj) (1 + o(l)), as n — 7- oo. Thus, by (g]) and for large n, Xj has normalized levels 



u 



nj 



UyAnr./T,) (1+0(1)). 



Consider u„(a;) = (itlj^^^^^-*-*, . . . , mIj^''''^'''''') = (nxi, . . . ,nxd) a vector of normalized levels of X and, for 
each j = 1, . . . , d, let {A"„j}„>i be an i.i.d. sequence with the same marginal distribution as {Xnj}n>i- 
Denote M„ = (Mni, ■ ■ ■ , Mnd) the vector of the componentwise maxima Mnj = V"=i j ~ ■ ■ ^ d. 

Proposition 2.4 We have Pxi in the domain of attraction ofGx., that is, 

lim P (M„ < u„(a;)) = Gx(a;i, Xd), 
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with 



= lim exp < 




/ n— >-cxD 





provided that Fy^ is in the domain of attraction of Gy{xi, ...jXd), where both Gx and Gy have unit 
Frechet marginals. 

Proof Just observe that 

d 

n 

Now the proof runs along the same Unes as in ([5]). □ 

In the foUowing, for any vector z and A C {1, . . .,d}, za denotes the sub- vector of z with indices in 

A. 

Example 1 // (Ti, . . . , Td) has support in {0, 1}'' and P (f]^^ j T, = 1, ^\^eD\.I Ti = = p{J), 7^ J C 
D = {1, . . . ,d}, then Gx corresponds to a geometric mean of the marginal distributions of Gy- In this 
case, we have 

Gx{xi,...,Xd) = Yl G'^'l\riXi,...,rdXd)j = W '^^■^ ( "TJT' ' ' ' ' "Tjf ) 

'p{{l})xi p{{d})xd 



n 

0#JCI5 



p{J) p{J) 



We illustrate the result with some choices for Gy ■ 

If the stationary sequence Y has common copula logistic, i.e. 



d 



Cy„ (iti, . . . = exp - I ^(-loguj) 

\3 = ^ 



then we obtain for Gy the logistic distribution and 



Gx(xi,...,xd)=exp|- (T.^f^J^^^)''^''] }' 

with 13 ji = p{{i}) / p{J) , which is an asymmetric logistic distribution already found in Tawn HSS}/ . 1990), 
by following a different probabilistic approach. The parameters /3ji increase the variability within the tail 
dependence coefficients regarding the departure distribution Gy- For X = {Xi, . . . ,Xd) with distribution 
Gx, we have, for instance, 

A^x.,x,)il,l)= limP(Fx.(X,)>l-l/t|Fx,(X,)>l-lA) = 2- ^ (pjl^'^+PJ-^" " 



where is the family of subsets of D containing {i,j}, expression that presents a larger number 

of possibilities than the value A(y. ^2 — 2°' of the symmetric logistic distribution. Note that 



Now, suppose that Gy{xi, . ■ . ,Xd) = H;^! Dfel-oo Aj=i cxp (-aa-ja;^^) , where {aikjj > 1, -oo < 
k < oo, I < j < d} are real non negative constants satisfying X]!=i Sfc=-oo '^ikj ^ I, j ~ 1, ■ ■ ■ ,d. This 
is the attractor MEV when Y is an M4 process (Smith and Weissman \32^ . 1996) and we have 



oo oo 



Gx{xi,...,xd)^Y[ n n A '^^p 

(=1 fc=-oo si^jcD jeJ 



and the bivariate tail dependence summarized by 



oo oo y- \ OO OO 



3 The multivariate extremal index 

The extremal index measures the tendency of clusters occurrence, a phenomena commonly observed in 
real data. In this section we will compute the multivariate extremal index of X„ (Nandagopalan |29j . 
1990). We start by analyzing some long range and local dependence conditions that will make easier its 
calculation. 

Proposition 3.1 // Y and T are strong-mixing, then X is strong-mixing. 

Proof Consider the events A and B, respectively, in cr-algebras generated by {Xi, . . . , Xp} and {Xp+s, Xp+g+i, . . .}, 
i.e., A £ i7(Xi, . . . , Xp) and B S cr(Xp-|_s, . . .). Given the independence between Y and T, we can state 

P{A nB) = E{P{A' n B'\Y)) = E (P(A'|Yi, . . . , Yp)P(S'|Yp+„ ...)) + 0{aT{s)), 

where A' £ cr(Yi, . . . , Yp), B' G a{Yp+s, ■ ■ •) and aT{s) denotes the mixture coefficients of the sequence 
T. The result follows by Doukhan [THI (1995, Theorem 3 in Section 1.2.2), since 

E {P{A'\Y,,. . . , Yp)P(i?'|Yp+„ ...))<E (P(A'|Yi, . . . , Yp)) E {P{B'\Yp+,, ...))+ 4aY{s), 

where ayis) denotes the mixture coefficients of the sequence Y. □ 

Consider u„(r) — •■• i "i^'''') a vector of normalized levels of X and Mi^„ =M„ = (M„i, . . . , M„d) 

the vector of the componentwise maxima M„j = Vr=i ' J ~ ■ ■ t As a consequence of the strong- 
mixing condition, we have 



lim P (m„ < u^-^)) = lim P f A/„i < u'-J^\ . . . , M„d < u^"^^^) 

n— ?-oc; \ / n— ^oo \ / 



[n/fc„] 



= exp <^ - lim nP Xi ^ u^"), f] X, < u^^) 



i=2 

where A:„ — > cxi and n/fc„ oo. as ri — > cxd. Therefore, X has multivariate extremal index 

- loglim„^oo P (m„ < ui'^^) 
dx{ri, ...,Td) = -rziz T—x-' 

- loglim„^oo P (M„ < n\r^ 

if and only if lim„^,oo nP [Xi ^ ul^\ M2.[„/fc„] < u^Zj exists and, if so, we have 



p(Xi ^ u1"\M: 



2,[ri/fc„] 



< U 



(t) 



6'x(Ti,...,Td) = lim — . (7) 

P (Xi ^ ui")) 
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This function relates the two multivariate extreme value distributions arising from the maxima in {X„}„>i 
and {X„}„>i. We analyze its values under additional assumptions on the scale sequence Y. 

Let 

rn I Tin \ 
and Uy^ (x) = {Uyh (2^1)7 • ■ • 7 {xd))- About the numerator in ([7]) we can write 



nP (Xi ^u(:),M2,[„/fc„] <uM) 




In some cases, the sequence Y is such that 

[n/k„ 



JimnP Yi^Uy. (^), f) Y, < Uy. ( J 



i=2 
k 



^limnP Y,^UY.(^),nY.<UY,(^) , 



4=2 

for some levels u[i'^ ~ Uy (7) and for some finite k, which makes easier the calculation of the multivariate 
extremal index. Situations leading to this easier procedure are highlighted in the next results. 

Proposition 3.2 Consider a = (ai, . . . , ad), b = (61, . . . , bd), Tnj with support in [aj, bj\, j = 1, . . . ,d, 



(^P^j >---7?^Yrf (^T^T^j j' /o*^ ^ ["j'^ji i = //Y satisfies condi- 

izonD^ (vi'^'''''\vi'^**'''^), defined by 

ln/k„] 

hm n V P ( Yi ^ v(-'(b)), Y, < v(-*(=^)), Y,+i ^ v'^'f^)') = 0, (9) 

i—k 

then D*'') ("""^)' defined by 

[n/fc„] 

hm n V P (Xi ^ u(7),X, < u(7),X,+i ^ uj^)) = 0, 

i — k 

holds for X, x = (ti, . . . ,Td). 

Proof Observe that, if we take ul^' = (t^Xi (n/Ti) , . . . , ?7xd (n/Td)), we have 

[«/fc„] 

hm n V P (Xi ^ uM,X, < u(:),X,+i ^ uf^: 

'i— /c 

= i^oy E ^ (yi ^ v(:'(^^)), Y, < v,(;*(^')), Y,+i ^ v(:'(^'+^))) dP(T„T„T.^,)(zi.2-z»+i) 

2— fc 
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Now, just observe that 



< p (yi ^ v^r Y, < vr*(=^», Y.+, ^ v(r(^») .□ 

Under the condition ([S]), smooth oscillations around v^i^ ''^^^ by Y (in the sense that its values which 
are at least lag-fc apart no longer exceed vl^ > vi^ ''^■''') are not followed by subsequent exceedances of 
vi^ (b))^ |.j^g Example SI we provide an illustration of such local behavior which leads to the condition 
D^'^) ('^"^^) ^'^^ scale model X. 

In order to avoid restrictions on the support of T^y , j = 1, . . . ,d, we propose now a greater restriction 
to the oscillations of Y around normalized levels. First observe that 



P (Yi ^ v(:'(^i», Y, < v^:* Y,+i ^ v(;'(^'+i))) < P (Yi ^ v(;*(^)), Y,+i ^ vi^' 



») 



with 



Yi 



Tin rdn 



^(zii V z,+i,i)'^iri ' ■ ■ ■ ' (zid V Zi+i,rf)/''iTrf^ 

Therefore, we can take for Y a sequence satisfying the following restriction on the local occurrence of 
exceedances of Vn which are lag- A: apart: 

[n/fc] 

limnY^P (Yi ^ vl-), Y,+i ^ vi^A ^ 0. (10) 

i—k 

This condition limits the size of the clusters of high levels and is satisfied by /c-dependent sequences Y, 
i.e., Y„ and Y„i are independent whenever |n — to| > k. 



Proposition 3.3 If Y satisfies ^^) for all v^T^ = Uyi (7), then X satisfies D^''^ {'^^^) ' f' 
quence of normalized levels u^^"* = Uxi (7). 



or any se- 



Corollary 3.4 //Y is k-dependent, then X satisfies D^*"'^ yui'^^j, for any sequence of normalized levels 

i-r) 

If X is strong- mixing and satisfies d''^' ('^""^) '^O'^'^ition, for all t G M^^, then it has multivariate 
extremal index given by 

P fXi < ul^\. . . , Xfc_i < uJr\ Xfc ^ u*"^ 



Ox{Ti,...,Td) ^ lim J — V , (11) 

P (Xi ^ u^") j 

whenever this limit exists (Chernick et al. [7] 1991, Ferreira [TT] 1994) . In this case, 

lim P (m„ < u^:') = cxp |- lim nP (x, < ui^\ . . . , Xfe_i < u^^), Xfe ^ u.^-^)) | . (12) 

The probability in the second term can be computed throughout the tail behavior of (Yi, . . . , Y^.). 
In order to improve such approach, define 

A(Y,,....Y.)(xi,...,Xfe)=^limtp('Yi ^Uy, ('^V...,Yfe^UY, (^))^ (13) 



with Xj G K^, j = 1, . . . , fc. This definition extends the upper tail dependence concept of Definition 2.3 
in Ferreira and Ferreira ([H], 2012). 

Note that, for with null components with j G {1, . . . , d}\Ji, 7^ 0, i = 1, . . . , fc, we have 

A(Yi,...,Yfc)(xi, . . . ,Xfc) = A((Yi).,^,...,(Yfc).,j((xi)ji, . . . , (Xfe),7j, 

where (xi)j. denotes the sub- vector of x.^ having components indexed in J,;. On the other hand, if x has 
some null component, we have Ayi(x) = 0. 

We are going to apply the result in (jl2p to our model and derive an expression for the multivariate 
extremal index, which will depend on the upper tail dependence function of (Yi, . . . , Y^) given in p3p 
and also on the dependence between random vectors Ti, . . . , T^. 

Proposition 3.5 // U3\) holds for all Xi, . . . , Xfc G ]R_^ then 



X Tf T T; 



k 



^0C/C{l,...,fc-1} \ / /U{A;}^ 



Proof Observe that 

nP fXi < u^;), . . . , Xfe_i < ui^), Xfc ^ 



n 



0C 

which leads to the result. □ 



05:i/C{l,...,/c-l} Vie/ / J 

/C{1,...,A;-1} \ie/U{/c} \^ / J 



Corollary 3.6 //X satisfies strong-mixing andD^ ' conditions, for all t G W^^, and if \1S\) holds 

for a// Xi , . . . , Xfc G then 

^x(n, . . . , r.) ^ 1 - V^-^^^^"-^-^^ ^ ^i^, (14) 

\0#/C{l,...,<i} \ / / , 



= 1 - > f E (-l)'^'"^A(^^,...,^^,,,^,, (75, . . . ) ) , (15) 



Moreover, for all x G Mf , we have 



Jim P(M„<Ux(nx))==exp<|-i?| (_i)m A(y,....,y. h^, • • • , 

3c/c{i....,fe-i} \ / /u{fc}^ 
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The result points out that, however much is the sequencial dependence in the stationary sequence T, 
the extremal index will be unit if Yi, . . . , are tail independent. This characteristic is illustrated in 
Example 21 with multivariate pRARMAX models, which are adjusted to a bivariate financial series in the 
last section. 



Example 2 Let {W„ = iW„i, . . . ,Wnd)}n>i be an i.i.d. sequence of random vectors with independent 
marginals Wnj, j = l...,d, having distribution function (d.f.) Fw^.(x) = (1 — and Ynj = 

Wn+ij V Wnj, J = f , . . . , d. Y is 2-dependent and thus we can apply the results on the calculation of the 

multivariate extremal index under the condition D^^^ (u^^M . We have A(Yi)j^(x/) — if \I\ > 1 and 



A(Y,.Y.)(xi,X2) = \im tP (y, ^ Uy, f^) , Y2 ^ Uy. (^^)) = J2 ^(^^ ^^2,) ■ (16) 



Then, by [I^ , 



^XlTl, ■ ■■,Td) 



r 



3 = 1 \ -J / 3 = 1 



1 - —E ( T!;"' A tP- 



Observe that the extremal indexes of the sequences {yry"}n>i, j ~ ^, ■ ■ ■ ,d, all equal to 1/2 while the 
presence of the random factors may increase this value for different 9j, j = l,...,d. Therefore, each 
marginal sequence may have a different tendency for clustering of high values. We also find that, since 
the marginals of Y„ are independent, the dependence structure of T„ does not affect the clustering of 
high values ofKn- expected from JT^P and illustrated in the examples of Section\^ only some tail 
dependence of Yn allows to account the dependence ofTnj, j = \, . . . ,d, on the value of 6x.{ti, . . . ,Td). 

Example 3 Suppose that Y is 2-dependent and {T„}„>i is a sequence of independent vectors with 
independent marginals having Bernoulli distribution with mean p. By we have 



^x(n, • . . ,Td) = 1 - - 



E p''^'(i-pri'^|(-iogGY.((^)^ 



0#JC{l,...,d} 

In the particular case of the previous example in U6]) . the function above becomes 

.2 d 



2 



9x{ri,. ■■,Td) 



d 

E-. 



and 6*^ = 1 - j = 1, . . . ,d. 

In the Example\^ ahead, we illustrate this choice ofTn to model lost values o/Y„, considering that 
the marginals Ynj, j ~ 1, . . . ,d, have total dependence. 
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4 Tail dependence 

The bivariate upper tail dependence of two random variables Xi and Xj can be measured through the 
tail dependence coefficient A(jCi,Xj)(l, 1) (Sibuya [21] 1960, Joe [TH| 1997). Extending this concept to 
cross-sectional lag-m, s upper tail dependence, m > 1 and < s < d, for a sequence X = {X„}„>i, we 
have 

X^^'''> = limt^oo tP > Ux, (i), Xi+,n,i+s > Ux,^^ (t)) . (17) 

If A^'*'' = we say that and Xi-f-m,i+s are upper tail independent. In this case, it is possible that 
a residual tail dependence captured at penultimate high levels may occur. This is measured through the 
asymptotic tail independent coefficient jy G (0, 1] (Ledford and Tawn [501 [H], 1996/97). Analogously, we 
can extend this concept to cross-sectional lag-m, s: 

P{Xi^l>UxAt/x),Xi+,n,l + s>Ux, + At/y)) X no\ 

hm — ^ — -, — - ~ hy '{x,y), (18) 

t^- P{X^,l>UxAt):X^+„^^+s>Ux,+At)) 

for all x,y >0, where ft,^'*'' is some non-degenerate function, homogeneous of order — 1/?7^'''\ 77^''''' G 
(0, 1] and such that ft,^"'*''(l, 1) = 1. In particular, we have 

P > UxAt):Xl+,n,l+s > Ux.^At)) = t-^/"-'"'''4"''^(0, (19) 

where ^^''^'^ is a slowly varying function, standing for the relative strength of dependence given a partic- 
ular value of r]'"x^'^K Whenever rf^'^^ = 1 and 4"'^-' converge to some constant < a < 1 then Xi^i and 
Xi+m,i+s are tail dependent (aI^'*'' = a), otherwise we have asymptotic tail independence with positive 
association if jy^"'"* > 1/2, negative association if r/^'*'' < 1/2 and (almost) independence if Ty^'*-* = 1/2 
(perfect if ^^'^^ = 1). 

By using the same arguments as in we obtain the following relation for \^^'^\ 
Proposition 4.1 The lag-m, s upper tail dependence coejficient ofX. is given by 

where K{Y-,^i,Yi+^,t+,){x,y) = YnTVt~^ootP {Yi.i > UY^[t/x) ,Yi+,n.i+s > Uy^+s (t/y)) is the lag-m, s upper 
tail copula of Y. 

Proposition 4.2 // (Yi_i,Yi^„i,i+s) satisfies il8\} for some function /i^"'"'' and coefficient riY^'^\ then 

p(xi,i > c/x,(t),xi+„,i+, > ux^^t)) = E ^h^^''''> (^^,^^±^^y-^/v^^^^^ 

for large t. 

Proof Just observe that, based on (??), and using ([T5|) . we have for large t 

P (Xi,i > C/xi (t) , > Ux^+s (t)) 

= P(yi,l > C7y,(t),yi+„,,i + , > Uy. + At)) / 4"''^ — 0'^m,i,Ti+„, + .)(^l,^2)(l + o(l)). 
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Now the result is straightforward from (jTH]) . □ 

Observe that the lag-m, s upper tail dependence or independence of X is ruled by Y. In the asymptotic 
tail independence case, this is even more evident since, by Proposition 14.21 they have the same lag-m, s 
asymptotic tail independent coefBcient, i.e., ?7^'*^ = 'r]'^''^\ only differing in the relative strength of 
dependence coming from the slowly varying functions. 

If we consider the models in Examples [5] and O we derive A^'*-* = 0, for all m > 1 and I < s < d. 
Observe that this corresponds to perfect independence of Y and thus it satisfies ^TE\\ and p^ . with 
rj'^''^ = 1/2, h'^''\x,y) = xy and l'^^^'\t) = 1. Therefore, 

P{Xis > Ux,{t), Xi+m,i+s > Ux,+At))^t-^ (21) 

i.e.. perfect independence of X too. 



5 Some particular cases that extend existing models 

In this section we consider our model restricted to a common rescaled factor Y, i.e., by taking Ynj = Yn' , 
7i > Oj J = I7 • • • J d, with Y^^^ = {Yn}n>i a stationary sequence having common Pareto-type d.f. with 
shape parameter a > 0, and r„j = j = 1, . . . , d, where Z = {(Z„i, . . . , Znii)}n>i is a stationary 
sequence, independent of Y^^^ with support M*^ and such that, rj = E{Z"j) < 00, j = 1, . . . ,d. We will 
therefore apply the previous results with /3j = a/'^j, j = 1, . . . ,d and AY(a;i, . • . , Xd) = Aj=i ^j- 

This case extends the heavy-tailed factor model in Lescourret and Robert ([22], 2006), which was used 
to measure the extreme dependence in loss severities of storm insurance data and where Y was considered 
a common latent factor corresponding to the intensity of the natural disaster. In addition, if we consider 
7 = 7j = 1, J = 1, . . . , d, we obtain an extended version of the scale mixture model of Li (US], 2009) with 
several applications to real data (see, e.g.. Arnold [3] 1983, Kotz et al. [12 2000 and references therein). 

As a consequence of Proposition 12.21 we obtain the following result for the upper tail copula of X: 

Ax(.„...,..)^i.| A^)=^(a^ 

for all {xi, Xd) 6 where {Zi, . . . , Zd) ==(Z„i, . . . , Znd)- From Proposition 12.41 we conclude that 
-Fxi is in the domain of attraction of 

G{x,,...,Xd)=c^p{~E\\/ --^L_I1L\ \ (22) 




Remark 5.1 // we assume that Z is an i.i.d. sequence of random vectors with i.i.d. marginals Znj, j 
1 . . . , d, having common d.f. H and rj ~ r, j ^ 1, . . . ,d, we easily derive 



and 



E {yU ^) = -r T.U ^ {im nti H ((../X,) V" z 



extending the bivariate results in Lescourret and Robert Q2^. 2006). 
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In the sequel, we rewrite Proposition 13. 2[ compute the extremal index and find the attractor MEV 
within this particular model. 

In what concerns the local dependence condition of Proposition 13.21 observe that 



\ l/a' 

n \ 



j=i \ y J , 



K>.v,-}^{v-,.A(t)"} 



and 

d da d 



Therefore, we obtain for this model the following particular result. 

Proposition 5.2 If Znj has support in [a, b], j — l,...,d, andY^^^ satisfies condition D^'^^ (^Vn^"^\vn'"''^^^ 
defined by 

ln/k„] 

lim 71 y p (Fi > < vi:^''^\Y,+i > i^-^w) ) = 0, 

i—k 

with r(a) = Vj=i ^ > Tib) = Vj^i > 0, j = 1, . . . , d, then dC^') (u^"^') holds for X, for each 

T = (ri,...,Td). 



Assuming that 



A(yi,...,n)(xi,...,a;fc) =^Iim tP (y, > f]^{UY. {t/x^)}^ 



(23) 



^ 

exists with (xi, . . . , x/e) E M^, then 



lim n P (Xi < ui^\ Xk-i < u^-^), Xfc ^ 



(t) 



3C/C{1,...,A:-1} = l ' 3 = 1 ^ 



with the convention that Ay^ ( Vj=i ) ~ Vj=i Therefore, if X satisfies strong-mixing and 

Hi^^ [vLn ^^ conditions, for aU r G M+, and if ^ holds for all (xi, . . . 6 1+, then 

(d ya d ya 

9x(ri, . . . , r,) = 1 ^ ^ ^ ^^^i^, (24) 
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the marginal extremal index is, for j — 1^ . . . ,d, 



J2 (-i)i^i+^v,.....nw>K.- 



/u{fc} 



(25) 



and, as n oo, P (M„ < Ux(?^x)) converges to the limiting MEV 



exp < -E 



z 



kj 



"-3 ' J 



7U{fe}/ J 



y0c/c{i,...,fc-i} 

for all (xi, . ..,Xd)£ K^. 

The cross-sectional lag-m, s upper tail dependent coefficient, to > 1 and < s < d, is given by 

^14 ^l+in,l + s 



(26) 



In case X is cross-sectional lag- to, s upper tail independent (i.e., A^'^"* = 0) and if (li, Yi+m) satisfies 
(fT5)l for some function Yi'™^ and asymptotic tail independent coefficient ry^T''', for large i, then 



^1,1 '^l+m.l + s 



t-^/"-"'4"^(t)(l+o(l)). (27) 



Next, we consider some particular examples. 



Example 4 [Multivariate autoregressive processes with random coefficients] Suppose that Y„ = c(y„_i V 
Wn), n > I, where < c < 1 is a constant and {W„}n>i is an i.i.d. sequence, independent of Yq, 
with common d.f F\y . The max-autoregressive sequence Y^^^ — {Yn\n>i, usually denoted ARMAX, has 
stationary distribution Fy if and only if < ~ Fw{x/c^)) < 00, for some x > and, in this 

case, Fy{x) = Fy{x/c)Fw{x/c) (Alpuim, ffj 1989). The Pareto distribution with parameter a > is a 
stationary df. 0/ Y with Fwix) = (1 — (ca;)^")/(l — x^"). 

If Znj has support in [a,b], j = 1, . . . , d, such that c < a/b, then Y'^^ satisfies condition D'^^ 

with T(a) and T{b) as in Provosition \5 .2\ Indeed, we have 



2/fe„ 



The first term is null since c < a/b 44> wi^'"'" < u^n'"^^^ / c and the second term is upper-bounded, succes- 
sively, by 

n ^P{Y, > v^^^^^PiW..^, > .(^C))) = n f (1 - F^.i^^^")) f 1 - -^^t^^^) 

Kn Kn \ FyA Vn / c) ) 

< ^{l-Fyivl^('^^)r^O 
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for any k„ 



oo, as n 



oo. Therefore, if the ARMAX sequence Y is such that c < a/b and Z has 



support in [a,b]'^, then X satisfies D^^^ [u„ ), for all t E M"? (Proposit 



From Ferreira and Canto e Castro (JT^, 2008; Proposition 3.6), we obtain the following bivariate 
upper tail copula function, for {x,y) G M^; 

A(yi,Fi+„)(a;,y) = a; Aye'"" 

Now we assume that Z is an i.i.d. sequence of random vectors with i.i.d. marginals Z„j, j — 1 . . . ,d, 
having common d.f H and rj ~ r, j ~ 1, . . . , d. After some simple calculations, we derive 



i.1 V^Ac-V^^^" 



i-YlH[{x,/x,y/'^z/ 

1=1 

d 



where r ~ rj, j ^ I, . . . ,d. 

Therefore, X has multivariate extremal index given by 



,Td) = l- 







+ E[^n{Z) [l-TTiZc)]) 


} 









where iriif) = Y^^i H {{tj /TiY^" , and for each j ~ 1, . . . ,d. 



E (Z" (1 - i?(Z/c)) + (cZ)" (1 - H{Zc))) 



In what concerns the limiting MEV ofX., for all {xi, . . . ,Xd) G , we have 



1 

lim P (M„ < Uxfrix)) = cxp { — V 



i=i 



E 



( 



\H{Z) 



TT* {Z)tT* {Z I C) 



{czr_ 

H{Z) 



Tr*{Z) [l-7r*(Zc)] 



where ^*{y) = Jlti H {{x./x^f/'^ y). 

Analogously, we obtain the cross-sectional lag-m, s upper tail dependent coefficient given by 



. (m,s) 



^-E {Z"H{c"'Z) + (c"'Z)" (1 - i?(c"Z))) 



Particular cases can be obtained by replacing H by any d.f. with support in [a,b], such that c < a/b. 

Now if we consider, Yn = Y^_i V Wn, n > I, with < c < 1 and {Wn}n>i on i.i.d. sequence with 
common d.f. Fw , independent ofYo, we obtain a max-autoregressive sequence with an exponent transfor- 
mation, denoted pARMAX (Ferreira and Canto e Castro, \12^ . 2008). A stationary distribution Fy exists 
if and only if < Ejli(l ~ Fwix'^')) < oo, for some x > 0, and is such that Fy{x) = Fy{x^^'^)F\y(x). 
This relation is satisfied by a Pareto(a) distribution with Fw{x) = (1 — .t^")/(1 — a;^"/"^) (see Ferreira 
and Canto e Castro JWj 2010 and references therein). The pARMAX sequence Y'^^ has lag-m upper tail 
dependence function K^Y^y^j^^^{x,y) = 0, for all m >1, and thus it is an upper tail independent process. 
Therefore, in this case, the lag-m, s upper tail dependence function ^(Xi i+s)ix,y) is also null, for 
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all m > 1 and s = l,...,d, the extremal index is unit and the limiting MEV is as in the i.i.d. case, 
i.e., it is given by Moreover, since Y^^-* satisfies U8\} with asymptotic tail independent coefficient 

r]^^ = 1/2 Vc™ and h^^\x,y) — a;yl{c'"<i/2} +2/"'^^^ !{£"»> 1/2}; x,y>Q (Ferreira and Canto e Castro, 
fW/ . 2008), by applying we obtain 



l{c-<l/2} + J7 



{c->l/2} 



' 1+s 



Observe that sequence X generated from an ARM AX recursion Y^^'' corresponds to a multivariate formu- 
lation of the RARMAX process introduced in Alpuim and Athayde (J^, 1990), with applications within 
reliability and various natural phenomena. If sequence X is generated from a pARMAX recursion Y, 
we have a multivariate formulation of the pRARMAX process introduced in Ferreira and Canto e Castro 
f\13f^. 2010), used in the modeling of financial series. 

Example 5 (Multivariate moving maxima processes with random coefficients) Consider Yn = 
Wn+i VWn), n> 1, where W = {Wn}n>i is an i.i.d. sequence with common d.f. Fw{x) = (1 — x^'^Y/'^ , 
X >\. The sequence Y'-"'^' is 2-dependent and thus X satisfies condition D^^''(u^), for all t £ M5|_, as well 
as, A(Yj Yj^^^) (x, y) = 0, for m>2. It is easily seen that 

\YuY^){x,y) = i (x + y-xVy) = ^ (xAy) 

Analogously to the ARMAX example above, we assume that Z is an i.i.d. sequence of random vectors 
with i.i.d. marginals Znj, j — l...,d, having common d.f. H and rj = r, j — \, . . . ,d, and a similar 
procedure lead us to 



E 



where r = rj, j — 1, . . . ,d. Hence, we derive the multivariate extremal index 
OxiTi,...,Td) = 





i-nti^((r.M)^/" z) 


) 




[T^UUH{{r,/r,y/'^Z)) 





the univariate extremal index, for j — 1, . . . , d, 

E{Z-[1-H{Z)]) 

C'Xj = 1 : 

and, for the limiting MEV, with {xi, . . . , Xd) G , 



Jim F(M„ < Ux(nx)) = exp j-i^x-^i; [ ^ \f{H ({x^ x,f"^ z)\ 



i=i 



The cross-sectional lag-1 , s upper tail dependent coefficient is 



1 
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The case m>2 corresponding to perfect independence ofY^^^ has been analyzed in \21\) . 

Observe that each marginal X n j is a moving maxima process with a random coefficient Z^ j, j — 
1, . . . , d. Therefore X corresponds to a multivariate formulation of this latter. For applications of multi- 
variate moving maxima processes, see Zhang (]S6^. 2009) and references therein. 

We are going to consider some particular cases of H . 



• (Multivariate Gumbel case) Suppose that H is a Frechet(5, S^) distribution, S,^ > 0, i.e., H{x) ~ 
exp{-Sx^^), X >0, with ^ > a. Therefore, P{Z°' < x) = exp(-(5x"«/") and EiZ") = e"/«r(l - 
a/S,) < CO, since ^/a > 1. We have 



1 



Y^ix./x.)^/'^ 



SY,{x,/x,)i/^ 



as well as 



E 



H{Z) 
Therefore, 



a 

\{h({x,/x,Y'^ Z 



r(i - a/0. 



G(a;i,. . . = exp |- ^^x, ^^"^ | 



which is a multivariate Gumbel distribution, and 



lim P(M„ < Ux(nx)) = exp |-2"/«-i (y^xJ^'"^ | • 

Thus the extremal index is constant and is given by 

ex(Tl,...,Trf)=2"/«-l=0„J = l...,d. 

• (Thinning model) Suppose that H is a Bernoulli(p) distribution, < p < 1. Since 

E ( ^ Tl R f (t: /t ^l/a Z n = ■n^^ _ r,^S^^=llf-i<-,} 



[H^^T{H{[x./x,Y/<^z)yp{l^p9 



and 



E 



H{Z) 



a 

X{h({x,/x,Y'^ Z 



p{l-pf^'='^<--^<-^\ 



G(xi,...,xrf)=exp<j -^x-i(l-p)^'-^f---.> 
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as well as 

lim P(M„ < Ux(nx)) =cxp<| - V ^-^(l - ^f^'<'=^' 
The extremal index is given by the function 



i=i 



?x(Ti,...,Td) = l- 



and 9j — I, j — 1, . . . ,d. 



6 An application to financial data 

The presented results besides give us an insight of the effect of the random factors on the extremal 
behavior of the model, suggest estimation methods whenever data for Y and T are available, through 
the estimation of the dependence in Y. Do not having real data sets in these conditions, we have chosen 
to illustrate in this section the application of the multivariate pRARMAX model considered in Example 
m to a bivariate stock market index data. The univariate pRARMAX was introduced in Ferreira and 
Canto e Castro ([13], 2010) and an adjustment algorithm to real data was implemented and applied to 
the S&P500 index, for the period April of 1957 to December 1987. Here we also consider the Dow Jones 
index for the same time period and our analysis is based on the volatility of the stock indexes measured 
through the square of the log-returns defined by Ri = log(fi+i)/ log(Pi), i = 1,...,7731 {n = 7731), 
where Pi is the i*^ closing day of the index. The values of i?,p and (i?^ ) , as well as, i?^'' and {R^'') , 
corresponding, respectively, to the log-returns and volatility of the S&P500 index and of the Dow Jones 
index, are plotted in Figure[TJ Observe that they have similar plots. The large peak corresponds to "Black 
Monday" stock market crash on the 19th October 1987. In Ferreira and Canto e Castro ([121, 2010) was 
considered a robust regression so that the marginals can be modeled through a Pareto distribution, and 
thus conducted the analysis on the transformed data X^^ ~ a (R^^) +b, with a = 13618.3 and 6 = 1.1. 

Following the same procedure, we consider X^'^ = a* (^R^'^Y + b* , with a* and b* estimated in 15154.2 
and 1.1, respectively. 

An adjustment algorithm for a pRARMAX process, Xn = YnZn, n > 1, obtained by considering 
r„ = y„^_i V Wn with < c < 1 and Fy{x) = 1 - x~^, where {Wn}n >i is an i.i.d. sequence independent 
of Yq and having distribution Fwix) = (1 — x~^)/{l — x~^^'^), and {Zn}n>i is a sequence of standard 
uniform random variables, was proposed in Ferreira and Canto e Castro f |13j 2010, Section 3). This was 
used in this latter reference to model the X^^ data producing the estimates c = 0.8 (although the values 
c = 0.85, 0.75 were also considered) and /3 = 2. Now we are going to apply this adjustment algorithm to 
the X^"^ data. In the following we summarize the running steps: 

1. Test if X^'^ is in the Frechet domain of attraction (see, e.g., Dietrich et al. [5] 2002) and estimate 
the tail index corresponding to /3~^ using, e.g., the Hill estimator (Hill, [17] 1975) and moments 
estimator (Dekkers et ai, [5] 1989). 

2. Estimate the parameter c, through the estimation of the tail index of = {S',^"''}, i — 1, . . . ,n ~ 1, 
where 



n + l-K^{XfJ)J \n + 1 - Kn{X^:{) J ' 
with Kn{x) = X)r=i ^{Xi<x}- 
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Figure 1: Left: log-returns of the S&P500 (top) and Dow Jones indexes (bottom); Right: volatility 
(i?f ) of the S&P500 (top) and Dow Jones indexes (bottom). 



3. Apply the criterion: "if Xf'-' > (X^^'^)'^, with the estimate c obtained in step 2, then X[''^ comes 
from the innovation component Wi " , to separate the innovations W and test if this sample is also 
in the Frechet domain of attraction. 

4. Capture the observations corresponding to Z, through the criterion: "if XP-' < {xPjiY and 
XP'' e ^ It: — ^ot-'^''^/w(t) < A ^ith fw the density function of W, 

TTi ^ P{Z, {XfJ^Y > W^)/P{{xPJ{y > Wi), TTo = 1 - TTi, and where /3 and c are the estimates 

obtained in steps 1. and 2., respectively, then, Zi ~ Xf''/ [X^'{Y "; v cannot be too large (larger 
errors) nor too small (not enough observations to carry out the test of the next step) . 

5. Test whether the sample of random variables Z captured in the step 4 has distribution Beta{j3 /c-\- 
1,1) (e.g., Kolmogorov-Smirnov test). 

For details, see Ferreira and Canto e Castro ([13] 2010) and references therein. 

In Figure [5] it is plotted the sample paths (against the k upper order statistics) of the heavy tailed 
test statistic (left), the Hill estimator (center) and the moments estimator (right), referred in step 1. 
The heavy tailed domain of attraction is not rejected for 50 < fc < 120 and 200 < fc < 1180, which are 
plausible values within this text (see again Dietrich et al. [S] 2002). In what concerns the sample paths 
of Hill and moments, the value where the paths yield approximately a flat line is at about 0.57 and thus 
we consider /? = 1/0.55 » 1.82. 

The Exponential and Pareto qq-plots and the empirical mean excess function plots (Figure [3]) indicate 
a Pareto- type model (see, e.g., Beirlant et al. [1] 2004). 

In Figure[5]we have the sample paths of the Hill (left) and the moments (right) estimators concerning 
the sample S (step 2). The plot is somewhat "stable" at about 0.8 and thus we take c = 0.8. 

Applying the criterium of step 3, we separate the innovations component W . The sample path of the 
heavy tailed test statistic (Figure O left) do not reject this assumption for k < 780, which is again a 
plausible result to conclude a Frechet domain of attraction. From the Hill and moments plots (Figure El 
center and right, respectively) we do not discard the same estimate of /3 « 1.82. 
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Figure 2: Left: sample path of the heavy tailed test statistic, against k upper order statistics, applied 
to the data X^^ (horizontal line: critical value above which reject the null hypothesis of a heavy tailed 
domain of attraction); Center and right, respectively: sample paths of Hill and moments estimators, 
against k upper order statistics, of the data X^^ concerning the estimation of the tail index /3~^. 
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Figure 3: Left to right, respectively: exponential and Pareto qq-plots of X^'' , empirical mean excess 
function against and against k upper order statistics. 




Figure 4: Left and right, respectively: sample paths of Hill and moments estimators, against k upper 
order statistics, of S concerning the estimation of c. 



According to step 4, we capture the observations corresponding to random coefficients, Z , considering 
the significance regions B^, for v = 0.05, 0.1, ... , 0.45 and c = 0.8. 

Implementing step 5, we applied the Kolmogorov-Smirnov test for the distribution Beta(l. 82/0. 8+1, 1) 
to the sample of random variables Z captured in the step 4. The rejection is obtained for v > 0.3. By 
Figure El the value v = 0.15 seems to be an appropriate choice (with 21 observations captured) and is 
in agreement with the simulation study conducted in Ferreira and Canto e Castro ([T3] 2010, Section 
4.1). Therefore, we consider that the model pRARMAX with c = 0.8 and f3 = 1.82, can be used for the 
modeling of the transformed data X^'' . 

Now we are going to estimate the dependence of {Xf^ , Xf-'), i — 1, . . . , 7731, based on the dependence 
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Figure 5: Left: sample path of the heavy tailed test statistic, against k upper order statistics, applied to 
the separated innovations W, according to step 3 (horizontal line: critical value above which reject the 
null hypothesis of a heavy tailed domain of attraction); Center and right, respectively: sample paths of 
Hill and moments estimators, against k upper order statistics, of W concerning the estimation of the tail 
index /3~^. 
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Figure 6: Empirical and theoretical distribution functions of the random coefficients Z captured on step 
4, for significance regions with, v = 0.05 ("u 05"),. . ., v = 0.45 ("u 45"). 



function A (Pickands [30] 1981) of a GEV distribution G, defined by 

G{xi,X2) = exp {-(.Tj'^ + X2'^)A{xi{xi + X2)^^)} ■ 

In order to have a dataset more coherent with an MEV model, we take block maxima within each 
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marginal. We use the Pickands ([30] 1981), the Capera et al. (0 1997) and the Hall and Tajvidi ( jl6j 
2004) nonparanietric estimators for the function A, and thus free of model assumptions underlying the 
data. We also consider the Lescourret and Robert estimator ([12] ■ 2006), given by 

max f(l - w)X,; ,wXi ^ 
A{'w) = min | l,max | — ^^^gp , w, 1 — w | | , (28) 



with 



/~SP ~DJ\ 
\Xi tXi \ 



-k:n 



-SP j 'I T DJ I I l{^.=^>X„-.;=:„=^X.°J>X„_fc,„°^}, 



where Xi:„-' < . . . < Xn-.n'^ , j G {SP,DJ}, are the respective order statistics. We have considered A{w) = 
{A{'w) + A{1 — w))/2 in order to constrain the estimators to be symmetric. This estimator was developed 
for factor models with a common factor Y , as described in Section [S] The multivariate pRARMAX model 
(Example 2]) was derived within this context. Observe in Figure [7] that the estimators are quite close, 
specially if we consider the bimonthly maxima, an indication that the multivariate pRARMAX may be 
a plausible model. 





Figure 7: Estimation of the dependence function A in ((28|) for the bivariate data [Xi^^ ^Xi^'^^, i = 
1, . . . , 7731, using the non parametric estimators of Caperaa et al. [5] (o), Pickands (+), Hall and Taj- 
vidi [H] (a) and using the estimator of Lescourret and Robert jH] (x; grey), considering (top-to-bottom 
and left-to-right, respectively) the weekly, monthly, bimonthly, three monthly, quarterly, semiannual and 
anual maxima. 
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